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We show that an effective magnetic-field acting on phonons naturaUy emerges in phonon dynamics 
of magnetic solids, giving rise to the phonon Hall effect. A general formula of the intrinsic phonon 
Hall conductivity is derived by using the corrected Kubo formula with the energy magnetization 
contribution incorporated properly. We thus establish a direct connection between the phonon Hall 
effect and the intrinsic phonon band structure, i.e., phonon Berry curvature and phonon dispersions. 
Based on the formalism, we predict that phonon could also have quantum Hall effect in certain 
topological phonon systems. In the low temperature regime, we predict that the phonon Hall 
conductivity is proportional to for ordinary phonon systems, while that for the topological 
phonon system has the linear T dependence with the quantized temperature coefficient. 
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I. INTRODUCTION 

The phonon Hall effect (PHE) has been discovered re- 
cently in Tb3Ga50i2(TGG)A!3: a magnetized solid could 
give rise to a temperature difference between two edges of 
the sample in the direction transverse to both the mag- 
netization and the driving thermal flow. The discovery of 
the Hall effect of the neutral carriers such as phonon has 
incited great theoretical interests^"—. Most of the theo- 
ries relate the effect to the Raman spin- lattice coupling^, 
and the standard linear response theory, i.e., the Kubo 
formula^^— or its equivalent*, are employed to calculate 
the thermal Hall coefficient. These investigations, while 
all focusing on the intrinsic limit, have not yet revealed 
a simple connection between the phonon Hall effect and 
the intrinsic phonon band structure for a general system, 
as that done in electron systems^. There is also a natu- 
ral and interesting question: could phonon systems have 
quantum (anomalous) Hall effect as well^ i^°i^^ ? 

In this paper, we establish a direct connection between 
the PHE and the intrinsic phonon band structure, i.e., 
phonon Berry curvature and phonon dispersions. To 
do this, first, we derive the general phonon dynamics 
applicable for magnetic solids, incorporating the Mead- 
Truhlar term in the Born-Oppenheimer approximation^^. 
The resulting dynamics contains an effective magnetic 
field acting on phonons, which gives rise to the PHE. 
It clarifies the microscopic origin of the spin-lattice cou- 
pling, and is readily amendable for the first principles 
calculation. Second, we calculate the thermal Hall coef- 
ficient of the system using the corrected Kubo formula, 
incorporating the contribution of the energy magnetiza- 
tio n^'^1^^ , which is overlooked in all previous calculations. 
As the result, we obtain a general formula for calculating 
the intrinsic phonon Hall coefficient. Based on the for- 
malism, we predict that phonon could also have quantum 
(anomalous) Hall effect in properly defined topological 
phonon systems. In the low temperature regime, we pre- 
dict that the phonon Hall conductivity is proportional to 
for the ordinary phonon systems, while that for the 



topological phonon systems has the linear T dependence 
with the quantized temperature coefficient. 

This paper is organized as follow. In Sec. II, we dis- 
cussed the general phonon dynamics in the magnetic 
field. In Sec. Ill, we presented our central formula for 
the phonon Hall conductivity. In Sec. IV, we gave the 
definition for the topological phonon systems and also 
discussed the possibility to realize it. In Sec. V, we show 
the law in the low temperature limit for the phonon 
Hall conductivity. In Sec. VI, we have a brief summary. 
We also include an Appendix section to show details of 
derivations. 



II. PHONON DYNAMICS OF MAGNETIC 
SOLIDS 

The Hamiltonian for the nucleus motion of a magnetic 
solid is determined by the complete form of the Born- 
Oppenheimer approximation^^: 



H = 



E 



i-ihVi^ - Ai^ {{R}))' 



^cff (R) , (1) 



where = d/dui^, and V,s {{R}) = {{R}) + 

Ei {{R}) is the effective interacting potential for nu- 
clei, including the direct Coulomb interaction between 
nuclei Ei {{R}), as well as the nucleus interaction 
mediated by electrons: Ec {{R}) — Eq {{R}) 

(nV2MK)EiK [(Vi.$o|Vi„$o) - |($o|Vi«$o) 
where Eq {{R}) is the energy of the ground state 
|<I>o({-R})) of the electron subsystem at the instanta- 
neous nucleus positions {R}. For the crystalline solid, 
the equilibrium positions of nuclei form a Bravais lattice, 
and we use {R} = + I = 1 . . . N; n — 1 . . .r^ 

with i?^^ = R^ + df^, where i?° is the center equilibrium 
position of Z-th unit cell, is the equilibrium position of 
K-th nucleus (with mass M„) related to the center, 
is the corresponding vibration displacement, N denotes 
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the total number of the unit cells, and r is the number 
of atoms in each unit cell. 

The most notable feature of Eq. (jlj is the 
presence of the vector potential ({-R}) = 

i/i($o({-R})|V,,$o({-R})) (V,, = d/dui^), as first 
pointed out by Mead and Truhlar— . In the modern 
language, the vector potential is related to the Berry 
phase^^. The corresponding "physical field" at the limit 
uii^ is: 



GS^'(i2°^i^?,)-2Mm 



5*0 



(2) 



where a, P = x, y, z, and the translational symmetry dic- 
tates that it must be a function of il° — iij! . 

In the external magnetic field, the ion will experience 
two vector potentialsi^: one from the real magnetic field 
and the other from the Berry phase. The role of the 
latter is crucial here. In a single atom, there is exact 
cancellation between the two, or the neutral atom would 
be deflected in the magnetic fleld. Therefore, it is not 
right to deal with the ion in the magnetic field directly 
by the "minimal substitution"^^. In the lattice system, 
the cancellation is not exact, and the spin-orbit coupling 
of the electron will give rise to nonzero GJJ^ — R'^r). 

To proceed, we adopt the usual harmonic approxima- 
tion by expanding Vcff ({-R}) to the second order of the 
vibration displacement ui^^^. For the periodic lattice, 
it is more convenient to use the Fourier-transformed dis- 
placement u^k = \/Mk/N ui^e'^^'^^'^ with k being 
the quasi-momentum, and the Hamiltonian can be writ- 
ten as: 
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k 



(PlPk + ulDkUk 



(3) 



where Pf, and Uk have 3r components with PK,a,k = 
-ihd / du,,a-k- A^a,k and Ut^a,k, a = a;, y, z respectively, 
A^aM = ih{^o\d^Q/du^a-k), and Dk is the 3r x 3r dy- 
namical matrix^*. 

The momentum P/^a.k has the commutation relation 
at the limit ^ 0: 



ihG 



kk' 



(4) 



where Gk is a 3r x 3r matrix with the component 

G'KC(,i//3(fc). 

The corresponding eigen-equation is: 



^kii^ki 



i 

-iDk iGk 



(8) 



where ipki is the ith eigen solution, and ujki is the cor- 
responding eigen frequency. For the non-hermitian Hk, 



we define ^ki = ipL^'k with Dk 



and 



Dk 

Isr X 3r 

-^3rx3r IS Unit matrix of dimension 3r x 3r, and ipki is 
normalized by ipkiipkj = Sij. Note that the 6r eigen- 
solutions can be divided into two groups of 3r positive 
and negative energy branches, and the two are related by 
cj^t) = and Vi;^ = ^/'itl* with 1 < « < 3r. This 

is a result of the symmetries = G_fc and = D^k- 
In the following we would only use the positive energy 
branches and drop the superscript for brevity. 

The Hamiltonian Eq. ^ can be diagonalized with the 
basis tpk- By defining the field operator ^k = (ufe, Pk)'^ , 
Eq. ^ can be rewritten as H = (l/2)X]fe*fc*fc with 
^k = ^k^k- Introducing the transformation 



3r 



*fe = E V f^kiipkiaki + \/haj-kiilJ*-ki^- 



ki ' 



(9) 



with 



a.ki,a 



kj 



^ijT hj — ^ ■ ■ ■ 3?", we can recover all the 



commutation relations, and diagonalize the Hamiltonian 
to: 



H 



3r 

E 

k;i=l 



hlUk 



ki I O^ki^k 



(10) 



Similar to the electronic dynamics in magnetic solids^-, 
the intrinsic phonon band structure is determined not 
only by the phonon dispersions, but also by the Berry 
connections of the phonon bands. We can define the 
phonon Berry connection as A-ki = ^ipki {dipki/ dk), and 
the corresponding Berry curvature as. 



Gna,K'p (fe) 



GS^'(i??)e-'^-(«?+''-') 



(5) 

where d^K' = — d^'. Gkq.k'^Cs) acts like an effective 
magnetic field for the phonon dynamics. Using Eq. ([2]) 
and ([5]), the quantity is readily calculable for the real 
materials using the first principles approach. 

The general phonon dynamics for a magnetic solid can 
now be determined. From Eq. (j3]), the linearized canon- 
ical equations of motion read: 



^ki = -Im 



"dk"" 



dk 



(11) 



We will show that the phonon Berry curvatures riki and 
the phonon dispersions ujki will fully determine the in- 
trinsic phonon Hall conductivity. On the other hand, the 
inter-band Berry curvatures proposed in some of previous 
studies are not needed in general^. 



III. PHONON HALL CONDUCTIVITY 



= , (6) Following the established general procedure^, we can 

; calculate the thermal Hall coefficient contributed by 

~ ~DkUk + GkPk , (7) phonons. For a magnetic system, the transport thermal 
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Hall coefRcient includes two parts of contributions: the 
usual linear response contribution k^'^^° and the contri- 
bution from the energy magnetization Me- 

(12) 



^.Kubo 



TV 



where Y is the total volume of the system, and T is the 
temperature. Me is the circulation of the phonon energy 
current, and the reason of the circulation can only be 
attributed to our effective magnetic field. 

A. Kubo contribution k^^"^" 



K^y^° is determined by the usual Kubo formula^: 



,Kubo 



hm^lini / dte-^*(j|_^; J|,j(t)) , 

(13) 

where JE,q is the Fourier transformed energy current op- 
erator in the wave-vector q, and (; ) denotes the Kubo 
canonical correlation^^. 

Using the procedure developed by Hardy^, and with 
the harmonic approximation and in the small q limit, we 
find that, [See Appendix Al]: 

JE^q ^l^^k{Vk + Vfe+q) ■^k+q + {h.c, q ^ -q) , 
k 

(14) 

where Vk = ^kHk, and (h.c, q — > —q) denotes the her- 
mitian conjugation of the first term after replacing q by 

With the phonon current operator JE,qy we have [See 
Appendix B[, 



where M-ki = Im 



k:i=l 



^ X 



and Uki = nsihhJki) 



dk ^ dk 

is the Bose-Einstein distribution. The presence of the 
"zero point" contribution (the extra ^ inside the paren- 
thesis) is due to the phonon number non-conserving 

terms (e.g., akiUk+qj or aL^fc+qj) energy current 

operator, which were often improperly dropped in many 
previous calculations^"—. As a result of the zero point 
contribution, K^y^° diverges when T ^ 0. 



B. Phonon energy magnetization 

The unphysical divergence can be removed by the sec- 
ond term of Eq. (|12p. The energy magnetization is de- 
termined as follows 

-'"'=Mi. (16, 



2M| - T- 



dT 



Ml 



1 



d 



h ■ 

'l — q, 'JE.q 



(17) 



where hq is the Fourier transformed energy density op- 
erator: 



K = \Y,^k^k+q' 



(18) 



The energy current operator defined in Eq. p4p should 
satisfy the scaling law necessary for the applicability^"^ of 
Eqs. (fT6HT7|l [See Appendix A2]. 

For the phonon energy magnetization in Eqs. (jl6| - |17p . 
we have [See Appendix C]: 



= ~2 ^ "^''^ 

k;i=l 



^ki^ki^'ki 



Ml, ( 2nkt + ujktn'k, + 1 



, (19) 



where n^^ — dnki/dujki- is obtained by integrating 
over the temperature T with the boundary condition that 
2M|; coincides with AI^ when T = 0. 



C. Phonon Hall conductivity 

We are on a position to present our central result for 
the phonon Hall conductivity. 



3h 



■ZphT ^ -f J dee^f^xy (e) 



dn (e) 
de 



where 



and 



(20) 



(21) 



3r 



^Ph - -j7 E ■ 

kA=l 



(22) 



Equation (|20p gives the general formula of the intrinsic 
phonon Hall conductivity for a magnetic solid. As ex- 
pected, the intrinsic phonon Hall conductivity is fully 
determined by the dispersions and the Berry curvatures. 

In the following, we show how the topological term 
emerges naturally in Eq. (|20p . 



First, we express K^y^°, Af|,'™*" and 7\,,f|^''"''^ as: 



^fy''"— I deemf: {e)n{e) , 



Mi 



M 



deem^^ (e) n (e) , 



(23) 
(24) 



z, intra ^ I j I 6r/\ 1 ~6r / \\ 2 (^) 

de [ m^, (e) - -ea^y (w) ) e 



(25) 
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where 



mil (e) =y (e - ^ki) , 



(26) 
(27) 



So the real energy magnetization is: 



Mh^-\j deemll (e) ^ d/3'/3' \2n (e) + e 



dn (e) 



1 



^ j dee'a% (e) ^ J rf/3'/3 



9e 



(28) 



where the boundary condition is that 2AI^ coincides with 
M|, when T = 0. There is no indefinity in Eq. (US]), for 
the thermodynamic quantity^ = Me/T should be 
zero in the high temperature limit. We can change the 
integration over 13' into the integration over e and we 
have: 

(29) 



1 

''2T 



^ TV ' 
dea% (e) 



dxxn [x) + e^n (e) 



(30) 



We can define a^y (e) = ^^^^fr-^, so we have: 



^xy 2x 



dxxn (x) + e^n (e) 



1 



, dn (e) 



The first term in Eq. ([3T|) is zero, for the following rea- 
sons: 



a% (e) = / dxa% ix) , 

J — oo 

<(-oo) = 0. 



(32) 

(33) 
(34) 

We can show Eq. ([55)) by the properties of the phonon 
Berry curvature: O^T-* — — fi'i^''^ where 1 < « < 3r, be- 
cause of the symmetry properties of ipki- Eq- P4p is zero 
because of the integration limit. Therefore, 



dn (e) 



I'y--^ / ^^''"'^ ^^^^ 



(35) 



Next, we express k*'^^ with the 3r positive energy bands. 



a% (e) = / dxa% (x) + I dxa% {x} , (36) 



Zph 
2^ 



(37) 



Using our definition in Eq. (|21l) . we can write Eq. ([57|) 
as: 



Therefore, 



dn (e) 



(39) 



Finally, after some simple algebra we obtain the 
phonon Hall conductivity in Eq. (PU)) and the topolog- 
ical term emerges naturally. 



IV. TOPOLOGICAL PHONON SYSTEM 

The first term in Eq. (PO)) is of the topological nature, 
determined by the global phonon band structures. This 
is different from the second term which is only deter- 
mined by the low energy sectors of the phonon bands 
limited by ksT. For the two-dimensional system (2D), 
is the Chern number, quantized as an integer. For 
the three-dimensional systems (3D), Zph is quantized— 
in the unit Gz/2tt, where Gz is the z-component of a 
reciprocal lattice lattice vector G. As the result, the 
phonon Hall conductivity has a topological contribution 
with the quantized linear temperature coefficient in the 
unit (7rfcs)V3/i (2D) or nklGz/Gnh (3D). 

Not surprisingly, most of the phonon systems have 
Zph = 0. It is natural to define the topological phonon 
systems as those with: 



Zph 7^ . 



(41) 



It puts a stringent constraint on what the real topolog- 
ical phonon system is: it requires that the sum of the 
Chern numbers of all phonon bands must be nonzero. 
The previous theoretical studies have successes in con- 
structing the phonon systems with non-zero Chern num- 
bers of individual band a^'^'-'i^^ . However, these models 
rely on the reorganization of the phonon bands within 
the positive energy branches, so still have zero Zph. They 
are not topological phonon systems in the stringent sense, 
nor will they manifest unusual phonon Hall conductivity. 
To realize the real topological phonon system, we need 
to look for insulating materials with spin-orbit coupling 
so strong that the resulting effective magnetic field can 
inter-mix and re-organize the phonon bands between the 
positive and negative energy branches. A realistic lattice 
model for realizing the topological phonon system is an 
interesting topic for future investigations. 

We note that the value of the topological contribu- 
tion to the thermal Hall conductivity is actually the 
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same as the longitudinal thermal conductivity of a dielec- 
tric quantum wire with Zph acoustic phonon modes^^— . 
This is because the topological phonon systems have chi- 
ral edge phonon modes that behave just like the ID 
acoustic phonon modes. The edge-bulk correspondence 
picture in the topological phonon system is very similar 
as that in quantum Hall systems. On the other hand, 
unlike the phonons in a quantum wire, phonons in the 
chiral edge modes can not be backscattered without cou- 
pling to the bulk. This will make topological phonon 
Hall effect more robust against imperfections than the 
quantized thermal conduction in ID quantum wire. The 
latter had been observed experimentally^^. 



LOW TEMPERATURE LIMIT 



where {r} = {ri,r2 ■ ■ - fN^} denotes the coordinates of 
the electron. We have: 

d<i>o {{R" + uo},{r}) 



duo 



■ exp 



J2 (p. - X n) $0 ({fi°}, {r - uo}) . (45) 



Substituting Eq. into (|43p . one can easily verify 
that the first term and the second term of the right hand 
side of Eq. exactly cancel each other due to the over- 
all charge neutrality. 



In the low temperature, only the low energy phonon 
modes are relevant for the ordinary systems. It is thus 
sufhcient to consider the long wave acoustic phonon 
modes. 



A. Constraint for the effective magnetic field 



An important constraint for this case is: 



^GS^'(i2°) = 0. 



(42) 



This can be verified directly by using the effective mag- 
netic filed in Eq. As we will show, this is still true 
even when we explicitly consider the external magnetic 
field because of the overall charge neutrality. It is easy 
to see that: 



d<fo {{R" + ttp}) 

duoa 



(43) 



where <I>o {{R^ + ua}) is the ground state wave-function 
of the electron subsystem when the whole system is dis- 
placed with Uq, and Ztot is the total charge number of 
the nuclei in an unit cell. We have assumed that there is 
an external magnetic field B. Since the displacement is 
equivalent to a re-definition of the origin, we have: 



$0 ({i2° + Mo},M) =exp 



-i- (B X no 



B. T'^ law 

It is easy to see that the phonon dynamics with the con- 
straint always has three acoustic modes that have zero 
energy at fc = 0, in consistent with the general require- 
ment of global translational symmetry. This remedies an 
important issue of the widely adopted phenomenological 
model of Raman spin-lattice coupling, which has non- 
vanishing coupling constant even in the long wave limit, 
inducing nonzero acoustic phonon energies at fc = 0. We 
will see that the constraint will change the theoretical ex- 
pectation of the low temperature behavior of the phonon 
Hall conductivity. 

With this constraint in mind, we can write down the 
general Hamiltonian for the long-wave acoustic phonons 
of an isotropic continuous medium^!: 



H = dx 



$o({fi°},{r-«o}) , (44) 



"(46) 

where P {x) — —\h5/5u{x) — ^[m], u{x) is the 
vibration displacement, /ii and /i2 are elastic con- 
stants, p is the mass density. The symmetry dictates: 
A[u\ = 71 W • (M X m) -I- 72 (M x u) , where M 
is the magnetization of the system, and 71 and 72 are 
coupling constants characterizing the spin-lattice cou- 
pling. With the Fourier transformation, we can identify 
Dap{k) = {1/ p){pi6ai3k'^ + P2kakf}) and Gap{k) = 
[-71 • M + (71 -f 272) fc^M^] with 
k — \k\ and a, (3,^ = x,y,z. We can then adopt 
our general formula for calculating the phonon Hall 
conductivity. 

To the first order of the magnetization M, we obtain 
the phonon dispersions: uj^i = c^k and U)f,2{3) — (^Tk ± 

72/cfe • M, where ct = \/pi/ p and cl = \/{pi + 2^2)/^ 
are transverse and longitudinal phonon velocities, respec- 
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tively. The phonon Berry curvatures are: 



51 (PM + kk-M) 

fc3 



(47) 



n 



fc2(3) 



k 92 {k^M + kk-M) 

fc3 P ' ^ ' 



where 51 - (71 + 272)(1 + 3J2)/(2ct^(52 - 1)), 52 ^ 
<7i5(3 + (52)/(2(l + 3(52)), j^j^d ^ ^ ^^/^^^ 

Using Eq. , we determine the phonon Hall conduc- 
tivity: 



_ 47r2fc|; 



71 + 272 + 5"^ +5 + 1 



(49) 

where we assume M is along the z-direction and the 
Debye energy hhj£, ^ /cbT. We can see that at the low 
temperature, k^^^ is proportional to T^, instead of T as 
proposed in the previous studies^. 

We can obtain insights from the above calculation on 
how the phonon would be deflected by the effective mag- 
netic field. We can see from Eqs. (|T7)) and pS)) that, 
on different branches, the phonon will experience differ- 
ent "reciprocal space magnetic fields". The correspond- 
ing anomalous velocity of the phonon is proportional to 
^ki X VT, similar to that for electrons^^. The net deflec- 
tion direction will be perpendicular to both the directions 
of the magnetization and the temperature gradient. 

A few comments concerning the disorder effect^^: In 
analogue to the anomalous Hall effect of electron sys- 
tems^, we generally expect the total phonon Hall coef- 
ficient can be decomposed to n^^y = k™ + n^^y + , 
where is the intrinsic phonon Hall conductivity we 



calculate in Eq. (|20|) , and the disorder will introduce the 
side jump contribution n^^y and the skew scattering con- 
tribution n^^™ ■ However, there is an important differ- 
ence between the phonon system and electron system: 
the mean free path of phonon in TGG had been de- 
termined to be ~lmm3'', much longer than its electron 
counterpart. Moreover, in the low temperature limit, 
the dominant contribution to the thermal conductivity 
is from the long wave phonons which don't "see" the dis- 
order^i. We thus expect that these disorder correction is 
less important in phonon systems, and the law of the 
phonon Hall conductivity will survive. 



VI. SUMMARY 



Hall conductivity of the ordinary phonon system is pro- 
portional to at the low temperature, while that for the 
topological phonon system has the linear T dependence 
with the quantized temperature coefficient. 



VII. APPENDIX 



A. Phonon energy current operator and its scaling 

law 



1. Phonon energy current operator 

We follow Hardy2^ to derive the phonon energy current 
operator. The Hamiltonian density is: 

/i(a;) = i^{A(a;-i2;,)i?,,+h.c.} , (Al) 



where = 1^ + V/„ and P/« = -i/iV/K - A^, 

and A(a;) is a localized function near a; = with 
/ da; A (a;) = 1. We adopt the Harmonic approxima- 

tion for Vu-. Vu = \ Ei'.'cp uuc^D^^p {Ri " R°) uv.'p 
where Z3^JD (i2° — is the dynamical matrix. The 

Hamiltonian is H — J dxh{x) = '^i^Hi^. The energy 
current operator is defined by the energy conservation 
equation: 



h{x) + V ■Je{x) = 0, 



(A2) 



where J e {x) is the phonon energy current operator, and 
h{x) = {l/\n) \h{x) , H . 



We can express 



h {x) , H 



as a divergence. We have: 



h{x),H\ = 4 E A(a;-i2i.) 



2h 

h.c. 



Ik.V K,' 



(A3) 



where " • • • " denotes the higher order terms such as 
' " Hifi, which will be a cubic product 



A (a; - Ri^) , 
of uik. and Pi'k,' . Noting that 



In summary, we establish the general phonon dynam- 
ics for the magnetic solids. Based on the dynamics, we 
propose a general theory of the PHE. Using the corrected 
Kubo formula, we link the intrinsic phonon Hall conduc- 
tivity to the phonon Berry curvature. The general for- 
mula suggests that phonon could also have quantum Hall 
effect, and our theory presents a rigorous definition of the 
topological phonon system. We predict that the phonon 



ih 



Ik 



2M«M^' 



(A4) 



and interchanging subscripts Ik and I'k' , we have: 
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= ifi [^{x- Ri.)- A{x- R, 



Ihi.l' . 



2M« 



Vi, 



(A5) 



Inserting the expansion 



A{x- Ri,) -A{x- Ri,.,) « - - Rt) 



dA{x~Rl,,) , dA{x-Rl 



dx 



dx 



into Eq. (IA5I) . we obtain: 

^\h{x),H ^V-Je (x) 



, (A6) 



(A7) 



where 



-A (0.-0) 



2M., 



, Vi> 



+ 



(A8) 



J 



Doing the Fourier transformation JE,q 
J dxJs (x) e~"''^, we have: 



JE,q = -^Ag {H^k [Dk + Dk+q) Uk^ 



-Pl^k [Gk + Gk + q) Pk + q , (A9) 



where Dk is the dynamic matrix with components 
the small q limit, Aq 1. We obtain: 



k 

+Vfcte + ^I+J)*fe+,. (AlO) 



Using the identities: ^kHk = Dk^kHk and ^kHl 
(Vfci?^) Dk, we obtain Eq. pi|. 



^. Scaling law for the energy current operator 
In order to calculate the energy magnetization, we need to verify that the energy current satisfies the scaling lawi^: 



jI (x) = (1 + {x)f Je (x) + O \v<'^^{x) 



in the presence of the gravitational field ■0 (r) that modifies the local Hamiltonian density by: 



(All) 



(x) = i^{(l + V^(a;))A(a;-i?,,)iJ,, + i/.c.} 



(A12) 



Similar to the derivation in Sec. I VII A 11 we have: 



V-jI{x)^^ 



k" {x) , H'f' 



(A13) 
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and: 



h"^ {x) , ijV-j = 1 / dx' (1 + V (x)) + {x')) r (a;, x') 



dx' 



dip [x) 



(1 + V {x)) r (a;, x') + (1 + ^ (a;)) {x' - x) ■ -h^T (a;, x') 

ox 



-i (1 + ^ (a:)) ^ (x;, - X,) K - Ir^r (a:, a:') + O [v^^'VX^) 



dxf^dxy 



with 



r(a;,a;')--^ ^{x-Ri.)[^{x' -Rvn'), 

Ik.V k' 



h.c. 



(A14) 
(A15) 

(A16) 
(A17) 



In Eq. (jA17p . we also ignore the terms which will lead to the cubic products of ui^ and i^/'^' in the energy current 
operator. For the first term in Eq. (jAlsp . after the integration over x' and repeating the derivation from Eq. (|A3[) to 
Eq. (|A7I) . we can show: 



- / da;' (1 + (a;))" T {x, x') = (1 + V {x)) V • Je {x) 



For the second term in Eq. (jAlSp : 



dx' (1 + V {x)) {x' - x) ■ ^t^T [x, x') 



i / , dtJj (x) 



dx 



A{x- Ri^) [r^,,, - Rl, [Hi,, H,^,] } , 



Ik.I' , 



.2{l + ^jix))^t^-JE{x). 



dx 

For the third term in Eq. (|Al6p . we have: 



I dx' (1 + ip (x)) Y K - « - X.) f^^r (a;, x') 



2h 



dxf^dxi, 
d'^ijj {x 



h.c. 



^(i + ^(a:)) j^^[^{^-Ri^)-^{^-Rv.')]{Rl.',-Rl,){Rl.'.-R'^^^^ 



dxndx- 



(A18) 

(A19) 
(A20) 

(A21) 

(A22) 

(A23) 

+ h.c. . 
(A24) 



Noting the expansion in Eq. Eq. (|A24D doesn't have ^^ere Vmj = 2 (tpkt^i'kj + V'L-^^Ij 

a contribution to the current at the long wave limit. 
Combining Eqs. (|A18HA2ip . we obtain Eq. (|A11|) . 



B. Derivation details for Kf^^° in Eq. ([T5]) 

Direct calculation of the Kubo formula leads to: 



■pki 
so: 



It is easy to verify 
dHk 



,Kubo 



V, 



E 



J — 1 



{iOki — ^kj) 



(A27) 



(A25) Note here i ^ so we have: 



9 



C. Derivation details for in Eq. (fT9)) 



— — — tlk- 



dk'j^ d ky dkj^ d ky 



dkx 



(A28) 



where we have used ipki^g^ - -^jr^vk], Vkj^^ 

-f ~ - 

~~aK"^ki^ and Hk = J2j ^kj^kjipkj and its hermi- 
tian conjugate. We can further simphfy the last two 

terms in Eq. ^M^ using = ^ - ^ipM and 



D 



OD 



k dky 



-gj^ — and we obtain: 



dkx 



HLDk 



dky 



-dk:'''''^ — 



dky 



dky 



dkj; dky dkj^ 

-^Hl^^k. 
dkj^ dky 



dkx 



Hk- 



dky 



ki ' 



The second term in Eq. (jA30l) 

,/,t dDufj dp^ji _ 
Vki dk^ "k Qf,^ Vki dk 

-'-^Hl^^M. and 
" 



hi dk^ 

Ji ~ 

dk^ 



dk^ 



kSK'^kt, and ^Hk-g^ 



k^k^k dk 

-1 



k dky ^k: 

fpki 



(A29) 



(A30) 



t 



for 

and then the two terms in Eq. (|A30[) 
vanish. Substituting Eq. (|A29p into Eq. (|A28p , we have: 



aofc ^ 

dka: 



dDj, 

dk^ 





{i^ki - ujkj) 

o T I ^'^ki fj dipkt 

=8wfc,;im I — — Hk- 



dkx 

Finally, we obtain: 



dk,, 



dipki fv dipkt 
-nk 



dk,, 



dkx 



(A31) 



6r 



.Kubo 



^ T ST- f '^'^ki fr '^'^ki dipki rr 9ll^k 

tm I — — tlk— — tlk 



2VT'" \ dkx dky 

k,i—l ^ ^ 
■ ^ki^ki , 

3r 



dky 



dkx 
(A32) 



ft -r / dipki rr ^Ipki dipki ~ dipk 

tm 1 — — tlk— —1 — tlk 



2VT 



k,i=l 

■ i^ki {2nki + 1 



d kx dky 



dky 



dkx 
(A33) 



From Eq. (jA32p to Eq. (jA33p , we have used the symmetry 
properties of ujki and ipki, so we come to Eq. (jlSp . 



After a direct calculation of the canonical correlation 
function in Eq. (jl7p . we obtain: 



ih d 



6r 



— 5^ ^ Sk+q,kjl {^k^k+qij + '^k+qMjl) 



k;i,j—l 



■UJkiLOk+qj 



^k-\-qj ^ki 



^ki — l^k+q 



, (A34) 



where Sk+q,kji = i^k+qji^ki and Vk.k+qij = 

V'fc.. / Ok V fc+qj+V^L ^ dk ^ k + qj- . 

First, we calculate the inter-band contribution from 
the terms with i ^ j m. Eq. (|A34p . When q — > 0, we 
have: 



Mi 



^ 6r^ Im (^^V'fejV^ij) WfciWfcj [Ukj - Uki 



8 ^ 

— 1 



^ki — <^kj 



(A35) 



We further have: 



Im 



V'fc^V^„ = -Im 



(A36) 

Noting V^jj = Vkji and inserting Eq. (IA36P into 
Eq. (jA35p . we obtain: 



6r 



16 



J — 1 



dky 



dky j '''^ 



OJki^kj {nkj — nki) 

^ki ~ l^kj 



, (A37) 



where we have interchanged i and j of the second term 
in Eq. (|A36| . Using Eq. (|A27| . we have 



32 



fc;i,j = l 



{uJki — i^kjY 



32 ^ 



fc;ij = l 



(A38) 



Comparing Eq. (jA38p and l|A25p . we finally obtain: 



3r 



^.^antcr ^ _ ( d'^kt jr ^'^kt _ dipkt ^ dipkt 



fc;i=l 

• Wfci (2nfci + 1) 

, 3r 



^ ^ ky d ky 



= -oYl ^l^'^k^ {2nk, + 1) 



dkx 
(A39) 

(A40) 



k:i=l 
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Second, we calculate the intra- band contribution from with lui\^^%hk,] = -Imf^/;i,,4^1 = lm]^^^hi. 
the term with i = j m hq. (|A34p : come to 



6r 



6r 



l^ki^k+qi (nk+qi — Tiki) 



OJki — ^k+qi 

For further simplification, note: 



q^O 



(A41) 



^, intra 
E 



fc;j:=l 



d {Hu - Ldkij 



dipki 



dkr 



dkx 



dky 



dk,. 



^L<^- (A44) 



dS, 



dqy 



^ {^k,k+qi + "^k+qM 



ki 



q^a 



dky 



ki ; 



dqy 



2 2V'fc: 



dkx dky 



q^O 



dkx dky dkxdky 



dkxdky 



(A42) 



nki 



l^k+qi ~ <^ki 



''ki 



q^a 



One can show that the sum of the last two terms in 
Eq. (|A42p is real, so it has no contribution. So we can 
write: 



6r 



M 



2, intra 
E 



-si: I- 



k;i=l 



„ d^Jki , dujki - dHk dtpki 
^^r^Vki— (- Vki 



dky 



dkx 



t 



dkx dky 



dkx d ky 



Using tpki- 

^ki 



(-Hfc-CJfci) 



9fex 



9±. 

"at. 



^ [Hk - ) and 
^ (^hI - ujkij , we obtain: 



6r 



M 



2, intra 



fe:i=l 



a-f/'fei f ~ \ d^ki 

Uk — ^ki 



dkr. 



d^ki 
dky 



(Hk 



- ^kj 



dij^ki 
dkx 



dky 



^fe»"L • (A45) 



Using the definitions for M.ki ^^"^ ^ki^ have: 



6r 



j^x, intra 
E 



k;i=l 
3r 



J2 i^l^+^k^^li)^l^^'k^ 



(A46) 

(A47) 



fc;i=l 



'kin'ki, (A43) FinaUy, combining Eqs. (|A40p and (|A47p . we come to 



Eq. ([19 
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